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Abstract Two alternative theories to dark matter are investigated by testing 
their ability to describe consistently the dynamics of the Milky Way. The first 
one refers to a modified gravity theory having a running gravitational constant 
and the second assumes that dark matter halos are constituted by a Bose-Einstein 
condensation. The parameters of each model as well as those characterizing the 
stellar subsystems of the Galaxy were estimated by htting the rotation curve of 
the Milky Way. Then, using these parameters, the vertical acceleration profile at 
the solar position was computed and compared with observations. The modified 
gravity theory overestimates the vertical acceleration derived from stellar kine¬ 
matics while predictions of the Bose-Einstein condensation halo model are barely 
consistent with observations. However, a dark matter halo based on a collisionless 
fluid satisfies our consistency test, being the best model able to describe equally 
well the rotation curve and the vertical acceleration of the Galaxy. 

Keywords Dark Matter ■ Milky Way ■ Modified Gravity • BEG Halos 


1 Introduction 

Astronomical observations on scales larger than several kiloparsecs indicate that 
the strength of gravitational forces cannot be explained only by the action of 
baryonic matter. A component, dubbed dark matter, with a cosmic abundance of 
about six times that of baryons and whose nature is still unknown, is necessary to 
explain different data as, for instance: i) the amplitude of peaks observed in the 
cosmic microwave background angular power spectrum; ii) the dynamical masses 
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of galaxy clusters, which are larger than those under the form of stars or/and hot 
gas; ill) the gravitational lensing effects produced by massive galaxies or clusters 
on more distant objects and iv) the formation of structures in the universe, the so- 
called “cosmic web”. However, in the literature, one of the most discussed effects 
is the “flat” rotation curve displayed by spiral galaxies that cannot be explained 
simply by the baryonic matter associated to the observed light distribution of 
those objects. 

Since the expected Standard Model relics do not have abundances sufficiently 
high to explain the aforementioned observations (excepting neutrinos, which have 
an adequate abundance but are not massive enough and have additional difficulties 
related to their free-streaming length), dark matter particles are expected to be is¬ 
sued from alternative models. Minimal Supersymmetric extensions of the Standard 
Model (MSSM) offer a plethora of candidates like gravitinos, photinos, s-neutrinos 
among others and presently, the neutralino, the lightest super-symmetric particle, 
is one of the most plausible candidates. So far, no signal of super-symmetry has 
been seen in experiments performed with the Large Hadron Collider [1]. More¬ 
over, results from direct search experiments are controversial. A positive signal 
modulated with a period of about one year is claimed to be present in data from 
experiments like DAMA/LIBRA, CoGent and CRESS-H [2l[3lll] that are not con- 
hrmed by the more sensitive experiments like XENONIOO and LUX [51161 . Negative 
results were also obtained from indirect searches related to the detection of high 
energy photons, leptons and hadrons originated from the annihilation of dark 
matter particles [ziEiiiiiniiiiiiii]. Thus, the only evidence to date for dark matter 
comes from its gravitational effects and, consequently, the possibility that General 
Relativity and the Newtonian theory break down at scales of kiloparsecs cannot 
be excluded a priori. 

This situation has stimulated the investigation of alternative modified grav¬ 
itational theories in order to explain, in particular, the “flat” rotation curve of 
galaxies. An early proposal in this sense, dating from more than 30 years ago is the 
Modified Newtonian Dynamics (MOND) [13]. This theory assumes the existence 
of a critical acceleration ao ~ 10~® cms~^. For motions in which the acceleration 
of bodies is much larger than ao, the gravitational force is essentially Newtonian, 
while in the opposite regime, the true acceleration is the geometric mean of the 
Newtonian value and oq. This simple theory is able to explain flat rotation curves 
and the Tully-Fisher relation m, as well as some dynamical aspects occurring in 
the local Universe m- However, MOND has difficulties to explain the dynamics 
of galaxy clusters as, for instance, the so-called Bullet-Cluster [T6| and the cosmic 
matter power spectrum Ha- 

Corrections to the Newtonian dynamics is a common feature of different models 
of effective low-energy quantum theory of the gravitational field HKIllSnj. Based 
on these ideas, the Renormalization Group approach was considered by |2ll|22l[23|, 
in which a ’’running” gravitational constant, depending on a given energy scale is 
expected. The consequences of the possible variation of G in galactic scales, under 
the assumption that the energy scale is fixed by the local Newtonian potential was 
investigated by the authors of reference |24| . According to them. Renormalization 
Group effects on General Relativity (RGGR) are able to explain the rotation curve 
of disk galaxies with a fit quality better than MOND or the Scalar Tensor Vector 
Gravity (STVG) [25]. Moreover, this modified gravitation theory hts quite well 
the observed velocity dispersion profile of elliptical galaxies |26| . 
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Another possibility often discussed in the literature assumes that dark halos 
are constituted by massive scalar fields or bosons that have undergone a Bose- 
Einstein condensation (BEC) in which they occupy the same quantum ground 
level Different investigations claim that the rotation curve of 

dwarf and low surface brightness (LSB) galaxies can be explained if a BEC halo 
is included in their mass distribution |291l33l[3l] . 

In the present paper the dynamics of the Galaxy is reviewed using either the 
RGGR theory or a dark halo constituted by a BEG. The Newtonian potential, nec¬ 
essary to fix the energy scale of the ’’running” gravitational constant in the RGGR 
model, was computed taking into account the contribution of different stellar sys¬ 
tems: the central bulge, the thin and the thick disks and also the neutral-|-ionized 
gas. The masses of the stellar components as well as Roo, a parameter charac¬ 
terizing the effective acceleration in the RGGR theory, were estimated by fitting 
the rotation curve of the Milky Way (MW). The same procedure was adopted for 
the BEC model. Once the parameters of these models were fixed by the fit of the 
rotation curve, the consistency of these theories was tested by computing and com¬ 
paring with data, the vertical acceleration at the solar position. Our calculations 
indicate that under these conditions the RGGR theory, despite the good descrip¬ 
tion of the galactic rotation curve, overestimates the vertical acceleration and this 
constitutes a potential problem for this model. The ht quality of the rotation curve 
resulting from the BEC model is worse than that of the RGGR theory but the 
predictions of the vertical acceleration are barely consistent with data, and alone, 
they cannot exclude this model. As we shall see, the difficulties for the BEC model 
concern the size of the halos as well as with the behavior of the rotation curve 
of the Galaxy beyond fS kpc. Only models of dark matter halos constituted by a 
collisionless fluid and having a mass distribution derived from cosmological simu¬ 
lations are able to satisfy our consistency test. This paper is organized as follows: 
in Section 2 the main aspects of the RGGR theory is briefly reviewed and the 
rotation curve of the Galaxy derived from this approach is discussed; in Section 3 
a similar analysis is performed for the BEG dark matter model; in Section 4, for a 
comparison with these models, the same analysis is made for a dark matter halo 
having a Navarro-Frenk-White (NEW) density prohle; in Section 5 the compari¬ 
son of all these models with vertical acceleration data is performed and, finally, in 
Section 6 the main conclusions of this investigation are given. 


2 The RGGR theory 

In this theory, quantum corrections produce a running gravitational constant G{fj.) 
that depends now on the energy scale fi of the theory (see EHil and references 
therein for details), i.e., 

= 1 + 2. Mm/mo) " 

where . is a small dimensionless parameter. We adopt here the same approach as 
|24| , namely, the energy scale is assumed to be related with the local Newtonian 
potential by a simple power law or, in other words 
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Using the field equations in the weak gravitational regime, it is possible to show 
that the effective potential of the theory in this approximation is 


4>ef 


4>n + 


2 Go 


(3) 


where 0jv is the usual Newtonian gravitational potential. In this case, the effective 
acceleration K^f, computed from the gradient of eql3] combined with eqslT]and[2l 
is given by 

d(j)N , d(j)N 


Since the expected circular velocity, assumed to be equal to the rotation velocity, 
is given by 


or 


= rK, 


ef,r 


(5) 


one obtains from the equation above and eqU] 






( 6 ) 


where Vat is the Newtonian circular velocity and we have defined = av(?. Note 

that the effective circular velocity depends on the Newtonian value corrected by 
a factor that depends on the local Newtonian potential and on the unique RGGR 
parameter Voo, which incorporates the product of the parameters, a and v. 


2.1 Application to the Milky Way 

The Galaxy is constituted by different sub-systems, which have their own mass 
distribution and kinematics. In order to evaluate the local Newtonian potential 
required to compute the effective circular velocity from eqlHl the following sub¬ 
systems were considered: the neutral-fionized gas, the bulge, the thin and the 
thick disks. 


2.1.1 The gas 


The projected density profile of the gas Sg{r) was taken from reference |35| . Since 
only the projected density profile is given, we assumed that the gas is distributed 
exponentially along the z-axis with a total scale of height 2Hg = 200pc. In this 
case, the gas density at the (cylindrical) coordinates r, z is given by 


Pg{r,z) 


Sg{r) /H„ 
2Hg 


(7) 


and the solution of the Poisson equation for an axisymmetric system, after inte¬ 
grating over the angular variable is 


(j)g{R,Z) = - 


4G 

Ha 


rSg{r)dr j dz 


K(u{r, z)) 


+r^ + 2rR+{Z - z)^ 




(8) 
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In the equation above K(x) is the elliptic integral of first kind and the function 
u(r, z) is given by 

AtR 

^ [i?2 + r2 + 2rR + {Z - z)^] 

Note that in the numerical computations, the values of the projected gas density 
given by m were multiplied by a factor 1.4 to take into account the presence 
of helium and trace elements. In this case, integrating the projected gas density 
given by [35] and correcting by such a factor, the total gas mass in the Galaxy is 
9.6 X 10® Mq. 


2.1.2 The bulge 


For the bulge, we have adopted the potential given by [36], i.e.. 


4>b{r,z) 


GMb 

(r2 + 2;2 + 52^1/2 


( 10 ) 


with b = 0.258 kpc. Fitting the very inner rotation curve of the MW and neglecting 
the possible dark matter contribution, the authors of reference [36] estimated the 
mass of the bulge to be Mb = 1.02 x lO^^M©. In our computations, we allowed 
the bulge mass to vary around this value when searching for the best ht of the 
rotation curve of the MW. 


2.1.3 The disk components 


The thin and the thick disks were supposed to have a double exponential mass 
distribution (radially and vertically) given by 


Pd{r,z) 


Md rIRi -zjHi 


( 11 ) 


In the equation above, Md, Rd and Hd are respectively the disk mass, the radial 
length scale and the scale of height of the considered component (thin or thick). 
The potential for a double exponential mass distribution was taken from reference 
m and, after some algebra is given by the expression 


4’d{R,Z) 


GMd dxJojx) [e - (^) e 




1 + ^ 

(1 ) 


( 12 ) 


where Jo(x) is the Bessel function of order zero. 

In our computations, the different length scales were kept hxed and taken from 
reference [S^, who have modeled the Galaxy using star counts in the J, H, and 
Ks bands derived from the 2MASS survey. These scales are respectively Rd = 
2.12 kpc, Hd = 0.205 kpc for the thin disk and Rd = 3.05 kpc and Hd = 0.64 kpc 
for the thick disk. 

In our fitting procedure of the rotation curve, the masses of both components 
were allowed to vary but not independently, since a constraint was imposed. In 
fact, the projected stellar mass density 17* at the solar neighborhood {tq = 8.3 
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kpc) is about 36 Mqpc~^ [39]. This value is essentially due to the contribution of 
both disk components (Note that including the gas contribution, the total baryonic 
projected density at the solar position is about 44 Mqpc~'^). Using ea llll the total 
projected stellar mass density is 

where the subscripts d\ and d2 refer respectively to the thin and thick disks. Thus, 
all the models satisfy the condition above. 


2.2 The rotation curve: RGGR model 


The circular velocity expected in RGGR theory was computed from eq|6| The total 
gravitational potential was evaluated by summing the contribution of the different 
components, that is 

4>N = + 06 + 4>dl + 4’d2 (14) 

whereas the Newtonian circular velocity was derived from 


V^ = R 


d4>N 

dR 


(15) 


In the fitting procedure, the masses of the bulge, thin and thick disks were 
allowed to vary, as well as the RGGR parameter Voo but not the gas distribution. 
Our best model is shown in figure 1 where the data points were taken from goj. 
These represent normal points issued from different tracers of the disk kinematics 
whose data were binned and averaged. 

This model is characterized by the following parameters: bulge mass = 0.78 x 
1O^°M0, thin disk mass = 0.98 x lO^^M©, thick disk mass = 2.60 x lO^^M© 
and RGGR parameter Voo = 226 kms~^. Hence, the resulting stellar mass of the 
Galaxy in this model is M* = 4.36 x 10 ^°Mq. 

The derived RGGR parameter corresponds to av = 5.67 x 10“^, which is about 
a factor 3.4 higher than that derived from the fit of the rotation curve of NGG 2403 
by the authors of reference jnj. They claim that the v parameter cannot vary from 
galaxy to galaxy but the a parameter can, contrary to MOND or STVG, which 
don’t have free parameters varying from one object to another. Despite the fact 
that the RGGR theory leads to a good fit quality of the rotation curve of the MW, 
the variation of the energy scale among galaxies is a weak point of this theory. 


3 The BEC model 

In this model, dark halos are constituted by an assembly of N identical self¬ 
interacting particles of mass m. In the Hartree approach, this self-gravitating 
structure is dictated by the uncorrelated single-particle stationary states of the 
mean-field potential created by the assembly of particles (Note that in the case of 
charged particles, the non-linearity of the eigenstates leads to a small correlation 
whose amplitude is of the order of the fine structure constant). The ground state 
corresponds to a condensed conhguration (BEC) in which all the particles occupy 
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Fig. 1 Best fit model of the rotation curve of the Galaxy using the RGGR theory. The total 
stellar disk mass in this model is 3.58 X 10^® the bulge mass is 0.78 X and the 

RGGR parameter is Voo = 226 kms~^. 


the lowest-lying orbital of the average potential. The Hartree approximation can 
be stated as the one-body self-consistent Schrodinger equation, in which the po¬ 
tential energy depends on the wave-function itself. If ^(r) is the single-particle 
wave-function and iZ'(r) is the system wave-function in the Hartree sense, then the 
particle density n(r) is 

n(r) =1 tf'(r) |^= N \ 4>{r) (16) 


and the single-particle wave function satishes the Schrodinger equation 

T) ^ 

- + mU{r)${r) = E^{r) 

where the mean held potential is dehned as 

U{r) = -GmN [ d^r' -h \ <P{r) 

Jbody \r-r' \ m2 


(17) 


(18) 


In the equation above, the hrst term on the right side represents the gravitational 
held and the second corresponds to the ground state of a ” hard-sphere” potential 
whose scattering length is a. 

When N ^ 00 , the solution of these equations can be obtained using the 
Thomas-Fermi (T-F) approximation, i.e., neglecting the so-called quantum cor¬ 
rection potential. In this case, the density prohle of the conhguration is simply 
given by 


P{r) = po 


sin{kr) 


kr 


(19) 
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The radius Rbec of the configuration is obtained with the condition kRbec = tt, 
corresponding to the position r = Rbec where the density goes to zero and is given 
by the relation 

The exact solution of Schrodinger equation, which can be obtained numerically, 
gives a density profile that goes only asymptotically to zero when r —^ oo (the 
wave function of the ground state has no nodes). In fact, with respect to the exact 
solution, Rbec includes about 92% of the total mass of the configuration. 

The gravitational potential due to the mass distribution given by eall9l is 


4>bec{r) 


GMbec r, , f Rbec \ . f 7rr \ 

Rbec [ ^ V J^^'^yRbec) 


( 21 ) 


Note that this equation is valid inside the BEC, i.e., for r < Rbec- When r > Rbec, 
the potential is simply given the point source approximation, this is 

_ GMbec 

0bec — 

r 

In these equations, Mbec = ^PoRbec/'’^ is the mass of the BEC in the T-F approx¬ 
imation. It should be emphasized that is important to distinguish the inner and 
outer solutions for the BEC gravitational potential because, as we shall see later, 
when fitting rotation curves the required BEC radius is often smaller that typical 
galactic dimensions. 


3.1 The rotation curve: BEC model 

In order to derive the circular velocity, the same procedure as before was used. 
Firstly, the total gravitational potential was computed, including all the MW com¬ 
ponents, i.e., 

(j}N = 4>bec + 4>g + (t>b + 4>dl + 4>d2 (23) 

where the same notation for the potential of the different galactic systems was used. 
Then, in a second step, the circular velocity in the galactic plane was computed 
from eall51 Again, the masses of the stellar components were allowed to vary in 
the htting procedure (observed the constraint mentioned before) as well as the 
BEC parameters. Note that a unique parameter appears in the RGGR theory but 
two parameters are required for the BEC model: the central density po (or the 
total mass) and the radius Rbec- 

The best fit of the rotation curve was obtained for a model in which the 
masses of the disk components are respectively 4.79 x 10^° Mq for the thin 
disk and 0.21 x 10^° Mq for the thick disk. The bulge mass in this model is 

1.1 X 10^° Mq. The BEC-halo has the following parameters: a central density 
po = 3.43 X 10^ MQkpc~^, a radius Rbec = 15.7 kpc and a mass of 1.69 x 10^^ Mq. 

The first point to be discussed is that in the optimization process of the rotation 
curve an additional constraint was necessary, that is the total disk mass was hxed 
with the value of 5.0 x 10^° Mq. This was necessary in order to have positive and 
non-zero values for the masses of both disk components. As a consequence, our best 
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Fig. 2 Best fit model of the rotation curve of the Galaxy using the BEG model. The total 
stellar disk mass in this model is 5.0 X 10^® the bulge mass is 1.1 X 10^® Mq and the halo 
is characterized by a mass of 1.69 X 10^^ Mq and a radius of 15.7 kpc. 


BEC model has a thin disk much more massive than the thick component. This 
is certainly a unrealistic result that could be avoided if the condition expressed 
by eall3lis relaxed. The second point refers to the radius of the BEC-halo, which 
is only 15.7 kpc, whereas galaxy-sized halos constituted by WIMPs are expected 
to extend up to 200-300 kpc. In fact, the analysis of the rotation curve of eight 
dwarf galaxies led to values for the BEC radius down to 1.0 kpc up to 12.6 kpc 
|33| . Although that investigation was focused on dwarf galaxies, it is interesting 
to recall that the BEC radius does not depend on the mass of the configuration 
(see ea l20ll but only on the scattering length and on the particle mass. Hence, we 
would expect that all BEC-halos should have similar dimensions but this is not the 
case. The third point concerns the derived mass of the BEC halo for the Galaxy. 
The mass of Galaxy (mostly in its dark halo) is still uncertain since values for 
its total mass included inside a radius of 100 - 200 kpc may vary by factors of 3 
to 4. In the one hand, using motions either of the Magellanic Clouds and of the 
Magellanic Stream [32] or the Sagittarius Stream [33] values around 5.5 x 10^^ Mq 
are obtained. On the other hand, the motion of the satellites of the Milky Way 
indicates higher values that are of the order of 2.0 x 10^^ Mq |441I451I36| . Thus, the 
derived mass of the BEC halo is about a factor of 3 less than the lower values and 
almost one order of magnitude less than the upper values of the estimated mass 
range. Finally, as it can be seen in hgure 2, the modeled circular velocity beyond 
15 kpc begins to decrease because the gravitational potential of the BEC-halo at 
these distances varies as 1/r as mentioned above. This is in disagreement with 
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observations, which seems to indicate a slightly decline of the rotation velocity 
only beyond 60 kpc m- All these aspects represent difficulties for this model. 


4 The NFW model 

In this section, the dark halo of the MW is modeled by a NFW profile, derived 
from fits of halo density profiles resulting from numerical simulations, in which dark 
matter is supposed to behave like a collisionless fluid. It should be emphasized that 
the computations presented in this section aim only to serve as a template to be 
compared with the previous models and to check our fitting procedure, since there 
is an extensive literature on this subject (see, for instance, |361I47| and references 
therein). 

The NFW density profile is given by the relation 


-1 



(24) 


In the equation above, the halo is defined by two parameters, pc and Rc, which 
are respectively a characteristic density and radius. The gravitational potential 
resulting from the solution of the Poisson equation with the mass density given by 
eg 1231 is 



(25) 


4.1 The rotation curve: NFW model 

The same procedure as before was adopted to compute the circular velocity, ex¬ 
cepting that in eg 12,31 the BEC potential was replaced by the NFW potential given 
by eol25l The best rotation curve was searched by varying, as previously, the mass 
of the stellar components and the two parameters defining the NFW halo. Fig¬ 
ure 3 shows the best fit model characterized by the following parameters: mass 
of the thin disk = 1.11 x 10^° Mq, mass of the thick disk = 2.52 x 10^° Mq 
and mass of the bulge = 1.07 x 10^° Mq. The parameters of the dark halo are 
respectively pc = 1.936 x 10^ MQkpc~^ and Rc = 17.46 kpc. These parameters are 
similar to those derived in reference m, whose authors performed a similar anal¬ 
ysis but including only the bulge and a single disk besides the halo and obtained 
Pc = 1.40 X 10^ MQkpc~^ and Rc = 16.10 kpc. The present parameters derived 
from the fitting of the rotation curve indicate that the halo mass inside 200 kpc 
is about 2.1 x 10^^ Mq, in agreement with upper bound estimates based on the 
dynamics of satellites of the Local Group |441l451l36| . 

The derived parameters of the MW halo implies that the expected dark matter 
density in the solar neighborhood is pdm{9>.2>kpc) = 0.75 GeVcm~^. This value is 
about twice the ’’canonical”’ value around 0.35 GeVcm~^ but is consistent with 
other recent independent estimates |36ll39ll3^ . 
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Fig. 3 Best fit model of the rotation curve of the Galaxy using the NFW model. The total 
stellar disk mass in this model is 3.63 x 10^^ the bulge mass is 1.07 x IO^^Mq and the 
NFW halo is characterized by the parameters pc = 1.936 x 10^ MQkpc~^ and Rc = 17.46 kpc. 



Fig. 4 Expected vertical acceleration from the considered models compared with data. 
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5 The vertical acceleration 

In this Section, the consistency of alternative models for dark matter, here rep¬ 
resented by a modified gravity theory (RGGR) and a dark matter particle model 
(BEG), is tested. Once the parameters characterizing the gravitational potentials, 
essentially the masses of the galactic subsystems, and those of the considered 
model are fixed by fitting the rotation curve of the Galaxy, they are expected to 
give also a correct description of the gravitational forces along the vertical axis. 
Thus, using the derived stellar masses for the galactic subsystems and the cor¬ 
responding parameters, the total potential can be estimated in any point of the 
Galaxy. Fixing the radial coordinate at the solar position, this is Rq = 8.3 fcpc, 
the vertical acceleration as a function of Z can be computed by the relation 

Kz = (26) 


5.1 The data 


Two sets of data were used in the comparison with expected values of the vertical 
acceleration shown in figure 4. The first is based on the literature, using reference 
|49| and on the sky surveys SEGUE and RAVE [51]. These are indicated by 
open circles in figure 4. The second set is based on computations performed using 
kinematic data on K-dwarfs taken from reference [39]. Their velocity dispersion 
data as a function of the distance to the galactic plane was fitted by a second order 
polynomial, i.e., 

o-z =-16.93-h 16.19Z+l.llZ^ (27) 


The fit is valid in the range 0.3 kpc < Z < l.l kpc and velocities are given in 
kms~^. The density profile n^{Z) of K-dwarfs along the vertical axis is given 
in references |391I52| and small corrections due to metallicity gradients along the 
vertical axis were neglected. The vertical acceleration, neglecting in the Jeans 
equation small terms coupling the radial and the vertical velocity dispersions, is 
given by m 


Kz 


1 d{nt,a%) 
n* dZ 


(28) 


Values of the vertical acceleration computed by such a procedure are shown as 
solid squares in figure 4. 

Simple inspection of figure 4 indicates that the RGGR model overestimates the 
vertical acceleration. The predicted acceleration, despite being able to reproduce 
adequately the effective force field along the galactic plane, fails to represent the 
gravitational acceleration along the the vertical axis. Concerning the BEG model, 
the predicted rotation curve gives a fit quality worse than that derived for the 
RGGR (or NEW) model but the predicted vertical acceleration agrees better with 
observation than the RGGR theory. 

Table 1 summarizes our results giving the score of these different models mea¬ 
sured by the reduced y-square, related either with the description of the rotation 
curve or the vertical acceleration profile. It is worth mentioning that the relative 
high values of the reduced associated to the analysis of the vertical acceler¬ 
ation are due only to one or two points whose observational errors were probably 
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Table 1 Reduced x-square for fits of the rotation curve and the vertical acceleration. The 
first column indicates the model, the second and the third give respectively the reduced 
and the degree of freedom relative to the fit quality of the rotation curve, while the last two 
columns give the same quantities relative to the vertical acceleration. 


Model 

xl/- 


xl/<^ 


RGGR 

1.83 

15 

- 

- 

BEG 

4.61 

13 

10.09 

14 

NFW 

2.17 

14 

5.27 

14 


underestimated. The RGGR theory gives a description of the rotation cnrve of the 
MW slightly better than the NFW prohle but overestimates largely the vertical 
acceleration. The BEG halo model gives the worse description of the rotation curve 
but represents the vertical acceleration better than the RGGR theory. Only the 
dark matter halo modeled by a collisionless fluid and a NFW density profile gives 
an adequate description of both the rotation curve and the vertical acceleration 
with the same set of parameters. 


6 Conclusions 

In this work, the consistency of two alternatives to dark matter was investigated. 
The hrst one concerns a modification of gravity, characterized by a running grav¬ 
itational constant (RGGR theory). The second one considers a model for galactic 
halos, which would be constituted by massive bosons forming a Bose-Einstein 
condensation. 

In both cases, the rotation curve of the Galaxy served as a departure point 
to £x the parameters of these two theories as well as those defining the different 
stellar sub-systems of the Milky Way like the bulge, the thick and the thin disks. 

Once the different parameters were fixed by the fit of the rotation curve, the 
consistency of these two models was tested by computing the expected vertical 
acceleration at the solar position. As it was shown, the predicted vertical accelera¬ 
tion profile by the RGGR theory overestimates the data derived from different sky 
surveys and this represents a potential problem for this model. Moreover, in the 
RGGR approach the variation of the ”a” parameter from one object to another 
seems to be necessary in order to explain the rotation curve of different galaxies, 
being an additional problem for this model. Our best BEG halo model gives a 
representation of the rotation curve of the MW worse than that obtained by the 
RGGR theory but the predicted vertical acceleration profile is barely consistent 
with data. However, the BEG model has also other difficulties. The radius of the 
BEG halo varies from galaxy to galaxy, contrary to the expectation of the theory 
since such a quantity does not depend on the mass of the considered object but 
only on the particle mass and on the scattering length. Moreover, in the case of the 
Galaxy, even taking into account the uncertainties in the present determination 
of its mass, the resulting value for the BEG halo is even smaller than the lower 
bounds. 

The usual dark matter model, in which halos are constituted by a collisionless 
fluid and whose density distribution is derived from numerical simulations, explains 
consistently both the rotation curve and the vertical acceleration in the solar 
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neighborhood. The halo mass inside 200 kpc derived from the fit of the rotation 
curve agrees with previous analyses and with other independent estimates based on 
the dynamics of the Local Group. The expected local dark matter density is about 
a factor two higher than the usual value of 0.3 GeVcm~^, but in agreement with 
other recent independent determinations. The present study indicates for the total 
baryonic mass of the Galaxy (bulge+disk+gas) the value of Mbar = 5.66 x 10^° Mq 
comparable with other recent estimates m- 
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